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Abstract. A theory of incompressible MHD turbulence recently developed by Boldyrev predicts 
the existence of a scale dependent angle of alignment between velocity and magnetic field fluctua- 
tions that is proportional to the lengthscale of the fluctuations to the power 1/4. In this study, plasma 
and magnetic field data from the Wind spacecraft are used to investigate the angle between velocity 
and magnetic field fluctuations in the solar wind as a function of the timescale of the fluctuations 
and to look for the power law scaling predicted by Boldyrev. Because errors in the velocity vector 
can create large errors in the angle measurements, particularly at small scales, the angle measure- 
ments are suspected to be unreliable except at the largest inertial range scales. For the data at large 
scales the observed power law exponents range from 0.25 to 0.34, which are somewhat larger than 
Boldyrev's prediction of 0.25. The results suggest that the angle may scale like a power law in the 
solar wind, at least at the largest inertial range scales, but the observed power law exponents appear 
to differ from Boldyrev's theory. 
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INTRODUCTION 

Phenomenological turbulence theories developed independently by Iroshnikov [Q]] 
and Kraichnan [|2|] extended and adapted the ideas of Kolmogorv's well known theory 
of hydrodynamic turbulence to incompressible MHD turbulence. Both Iroshnikov and 
Kraichnan predicted an equipartition of energy between kinetic and magnetic field 
fluctuations in the inertial range and an energy spectrum proportional to k~ 3//2 . But 
these early studies neglected the anisotropy of the turbulence which was subsequently 
found to be ubiquitous in both laboratory plasma experiments and in theoretical studies 
based on analysis and simulations of the equations of resistive incompressible MHD dH]. 
Turbulence in magnetized plasmas is spatially anisotropic with the local mean magnetic 
field providing a natural axis of symmetry. 

More recent theories of incompressible MHD turbulence incorporate the anisotropy 
of the turbulence into the theory in a fundamental way. An influential theory of this kind 
is the theory of Goldreich and Sridhar (1995) [4], hereafter GS95, with important correc- 
tions by Goldreich and Sridhar (1997) [5]. GS95 introduced the idea of 'critical balance' 
in which there is a balance between the eddy turnover time, or energy cascade time, and 



the Alfven crossing time of two wavepackets propagating in opposite directions along 
the local mean magnetic field. As a consequence of critical balance, the GS95 theory 
predicts that for a typical wavepacket the wavelengths parallel and perpendicular to the 

local field satisfy the anisotropic relation Aj_ Am and the perpendicular energy spec- 
trum is proportional to k ± 5 ^ 3 . 

The decade following the publication of GS95 saw improved simulations of incom- 
pressible MHD turbulence in two and three dimensions, many of which showed that 
for plasmas with a strong mean magnetic field |Bq| comparable to or greater than the 
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r.m.s. magnetic field fluctuations the perpendicular energy spectrum scales like k ± 
in contradiction to the GS95 theory la III, Is |9fl. To resolve this discrepancy Boldyrev 
EEl developed a new phenomenological theory with a critical balance condition dif- 
ferent from that of GS95. A somewhat different theory was derived by Beresnyak and 
Lazarian [12]. The key new idea introduced by Boldyrev is that as energy cascades from 
large to small scales through the inertial range the velocity and magnetic field fluctua- 
tions undergo an alignment process whereby the average angle between 8v± and 8b± 
is a monotonically decreasing function of scale. In his theory, Boldyrev predicts that 

the angle obeys the scaling law °c Aj^ 4 and that the perpendicular energy spectrum is 
proportional to k^ 1 . 
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Evidence for Boldyrev's alignment process and for the scaling law 6 A ± have 
been obtained through direct numerical simulations of forced, steady state incompress- 
ible MHD turbulence in three dimensions by Mason et al. Il3l [Till. These simulations 
demonstrate that for a number of different types of forcing functions the system develops 
an inertial range spanning approximately one decade in wavenumber where the perpen- 
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dicular energy spectrum is proportional to k ± and, presumably over the same range, 

~ Aj /4 . The successful demonstration of Boldyrev's alignment process by Mason et al. 
prompted us to consider whether this alignment process may also take place in the solar 
wind, a naturally occuring turbulent plasma that is directly accessible to in-situ space- 
craft measurements [ 3 [la [l^ Oil]. Therefore a study was undertaken to investigate the 
possible existence of a scale dependent alignment between velocity and magnetic field 
fluctuations in the solar wind. The results of this study shall now be briefly summarized. 



THE QUANTITY BEING MEASURED 

Boldyrev and his colleagues measure the angle using the formula lfl3l.[I 



where 8v± and 8b± are the projections of the fluctuations 8v = v(x + r) — v(x) and 
8b = b{x + r) — b(x) onto the plane perpendicular to the local mean magnetic field 
Bq(x), respectively, and the displacement r is perpendicular to Bq(x). It is desirable 
but not possible to employ precisely the same formula (1) when analyzing solar wind 
data. A single spacecraft essentially performs measurements along the solar wind flow 



direction. Because the flow is super-Alfvenic, Taylor's "frozen turbulence" hypothesis 
may be used to relate the time T between measurements to a spatial separation r = V SW T 
along the average flow direction, approximately the radial direction in heliocentric 
coordinates. The angle between the mean magnetic field and the average flow direction 
varies significantly due to variations in the direction of the mean magnetic field, but 
it is often near 45 degrees, the inclination of the Parker spiral at 1 AU. Even though 
the theory is formulated for fluctuations between two points x and x + r inclined at 90 
degrees to the local mean magnetic field, we expect that any alignment (if it exists) will 
still be measureable at inclinations of 45 degrees lfl9ll . 

The solar wind data employed in this study consists of simultaneous measurements 
of the average velocity vector and magnetic field vector from instruments on the Wind 
spacecraft. The data was acquired near the orbit of the earth at 1 AU when Wind was far 
away from the influences of the Earth's magnetosphere and bow shock. The data have 
an approximate 24-second cadence. The intervals analyzed have durations from several 
days to several months. For a given time scale T the fluctuations 8v = v(t + t) — v(t) 
and 8b = b(t + t) — b(t) are projected onto the plane perpendicular to the local mean 
magnetic field Bo(t) to obtain 8v± and 8b±, respectively. The average angle d(r) 
is then computed using equation (1). The local mean magnetic field Bo(t) is scale 
dependent and is defined as the average (vector) over the interval from t — x to t + 2x. 
The essential difference between this study and the method employed by Boldyrev is 
that the displacement r is not perpendicular to Bq in the solar wind. 



RESULTS 

An example of the angle measurements obtained using Wind data is shown in Figure 
1 . Note that the measurements span most of the inertial range, the range of timescales 
from approximately 3 s to 10 4 s. All the data analyzed in this study show an approximate 
power law scaling at the largest inertial range scales (a straight line on a log-log plot) 
followed by a roll-over of the curve at medium to small scales. This roll-over is clearly 
a deviation from the power law behavior seen at the largest scales. We suspect that the 
roll-over is due to small errors in the plasma velocity measurements which can cause 
large errors in the angle measurements at small scales. For this reason, estimates of the 
power law exponent p in the scaling law G(x) %p are only performed at the largest 
inertial range scales from 2 x 10 3 to 2 x 10 4 s. A linear least-squares fit on a log-log plot 
yields the best fit line and power law exponent p = 0.34 for the data shown in Figure 1. 

Proton velocity moments obtained using the 3DP instrument on Wind are subject to 
measurement errors, like any measurement apparatus. Let e denote the typical measure- 
ment error for the velocity components v x , v y , and v z . Then, the typical error for the 
velocity difference 8v x is 2e. Suppose the errors in the magnetic field 8b are negligi- 
ble compared to the errors in the velocity 8v. To analyze the error, choose the z-axis 
parallel to the vector 8b and let the vector 8v lie in the jcz-plane. Then an error 2e in 
the measured velocity 8v x will create an error 80 in the angle 9 as determined by the 
equations 

rin(e)=** sin (e + 50) = ^^. (2) 
ov ov 
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FIGURE 1. The average angle 9 as a function of the timescale T (circles) obtained from Wind 3DP 
data. The solid red line is the best fit line over the range from 2 x 10 3 to 2 x 10 4 s. The slope of the line 
yields the power law exponent p = 0.34. Each data point represents an average over many measurements 
(~4xl0 5 ). 



Assuming the angles are all small, it follows that 89 ~ 2e/|5v|. If |5v| is replaced by 
its r.m.s. amplitude 8v, then the quantity 89 = 2e/8v can be thought of as the smallest 
resolvable change in angle. 

This rough error estimate can be used to derive rough error bars for the angle mea- 
surements in Figure 1 . We emphasize that these are not true error bars, only a rough indi- 
cation of the smallest resolvable angle at each scale. Assuming a typical value e = 0.76 
km/s, the rms velocity fluctuations 5v(t) are computed from the data and used to com- 
pute the estimated error 89 = 2e/8v. The result for 89 is shown by the dashed orange 
line in Figure 2. Note that at the smallest scales the magnitude of 89 is comparable 
to the measured average angle 9. This leads us to believe that the measurements at the 
smallest scales are unreliable. The error bars in Figure 2 furthermore suggest that data 
in the middle of the range are unreliable too. Only the data at the largest inertial range 
scales are considered to be sufficiently reliable for the purpose of estimating the power 
law exponent in the scaling law. 

Results for the power law exponents estimated from linear least-squares fits over the 
range of timescales from 2 x 10 3 to 2 x 10 4 s are shown in Table 1. On average, the 
values for the power law exponent p obtained from the data are larger than the value 
0.25 predicted by Boldyrev's theory. Therefore, the measurements are not in overall 
agreement with the theory. However, in two of four cases studied, intervals 2 and 3, the 
values p = 0.25 and p = 0.27 are within 8% of the predicted value 0.25. This prompts us 
to ask whether there may be differences in the solar wind conditions for intervals 2 and 3 
compared to intervals 1 and 4, differences that may make it more likely to find agreement 
with the theory in the former case than in the latter. One difference is that intervals 2 and 




10 



10 10 10 
Timescale (s) 



FIGURE 2. The average angle 9 (circles), the error 89 — 2e/8v assuming e = 0.76 km/s (dashed), and 
the corresponding error bars 9 ±89 for the same data shown in Figure 1 . Note that these are not true error 
bars, only a rough indication of the smallest resolvable angle at each scale. 



3 consist of nearly homogeneous low-speed wind, devoid of any high speed streams, 
while intervals 1 and 4 both contain a notable presence of high speed streams. Recent 
work by Dasso et al. [20] has suggested that low-speed wind appears to be dominated 
by wavevectors perpendicular to the local mean magnetic field, as expected for quasi-2D 
turbulence, whereas high-speed wind appears to be dominated by wavevectors parallel 
to the mean field. If the turbulence in intervals 1 and 4 is not predominantly quasi-2D 
turbulence, then Boldyrev's scaling law would not be expected to apply and that may 
explain why we observe deviations from Boldyrev's scaling law in these cases. This 
interesting possibility needs further study. 



CONCLUSIONS 



The study of solar wind turbulence presented here suggests that the angle appears to 
scale like a power law at the largest inertial range scales, but in many cases the observed 



TABLE 1. The power law exponent p obtained from fits of 
the angle 9 over the range from 2 x 10 3 to 2 x 10 4 s. 



Interval Begin Date 



End Date 



Days 



P 



2 
3 
4 



01 Jan 1995 29 Jul 1995 

15 May 1996 16 Aug 1996 

08 Jan 1997 09 Jun 1997 

23 Aug 2000 15 Feb 2001 



209 
93 
152 
176 



0.34 
0.25 
0.27 
0.34 



power law exponent is significantly larger than the value 1/4 predicted by Boldyrev's 
theory. There is a clear breakdown of the scaling behavior at medium- to small-scales. 
Whether this breakdown is partly a genuine physical effect or is solely caused by 
measurement errors in the velocity data is not known at the present time. Therefore, we 
do not know whether the apparent power law scaling seen at large inertial range scales 
continues to hold at medium- and small-scales. More precise plasma measurements with 
errors less than 0. 1 or 0.05 km/s for each velocity component are needed to answer this 
important question. 

In half of the cases studied the power law exponent was within 8% of the value 0.25 
predicted by Boldyrev. In these cases the solar wind can be characterized as generally 
low-speed wind without any high-speed streams. The better agreement with Boldyrev's 
theory in these cases may be because low-speed wind is more quasi-2D in nature than 
high-speed wind as suggested by the analysis of Dasso et al. [|20|]. This is an interesting 
open question. 
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